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ABSTRACT 

An improved convection heat-transfer model has 
been developed for the prediction of the transient 
tooth surface temperature of spur gears. The dissipa- 
tive quality of the lubricating fluid is shown to be 
limited to the capacity extent of the thermal boundary 
layer. This phenomenon can be of significance in the 
determination of the thermal limit of gears accelerat- 
ing to the point where gear scoring occurs. Steady- 
state temperature prediction is improved considerably 
through the use of a variable integration time step 
that substantially reduces computer time. Computer- 
generated plots of temperature contours enable the 
user to animate the propagation of the thermal wave as 
the gears come into and out of contact, thus contribut- 
ing to better understanding of this complex problem. 
This model has a much better capability at predicting 
gear-tooth temperatures than previous models. 

NOMENCLATURE 

A heat flux area, m 2 (in. 2) 

a acceleration of fluid particle, m/sec 2 

( in . / sec^) 

b Herzian contact width, m (in.) 

C heat capacitance matrix 

C a dimensionless constant (2R 0 /F) 2 , see Eq. (9) 

F face width of tooth in contact, m (in.) 

f coefficient of friction 

G centrifugal acceleration, m/sec 2 (in. /sec 2 ) 

hj heat-transfer coefficient for lubricated flank 

of gear, W/hr m 2 K (Btu/hr ft 2 °F) 

hj t heat-transfer coefficient for top land of 
gear, W/hr m 2 K (Btu/hr ft 2 °F) 


h s heat-transfer coefficient for gear sides, 

W/hr m 2 K (Btu/hr ft 2 °F) 

h t heat-transfer coefficient for unlubricated 

flank of gear, W/hr m 2 K (Btu/hr ft 2 °F) 

J heat-conversion factor 

k thermal conductivity, W/m K (Btu/ft °F) 

N number of teeth 

n normal coordinate to tooth profile, m (in.) 

diameter to pitch ratio, m - ^ (in.”l) 
q heat flux, W/hr (Btu/hr) 

R 0 addendum radius of gear, m (in.) 

r radius vector of fluid particle, m (in.) 

S surface finish 

s tangential coordinate to the tooth profile, m 

(in. ) 

T temperature, K (°F) 

t time, hr 

V s instantaneous sliding velocity, m/sec (ft/sec) 

V t instantaneous total velocity, m/sec (ft/sec) 

W normal tooth load intensity, N/m (lb/ft) 

x,y,z dimensionless Cartesian coordinates of gear 
tooth 

a thermal diffusivity, m 2 /s (in. 2 /s) 

y finite difference weighting factor (0 < y < 1) 
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0 dimensionless gear temperature 

p dynamic viscosity 

v kinematic viscosity 

w angular velocity vector, sec”* 


Considering a frame of reference, attached to the 
gear at its center of rotation and thus rotating with 
an angular velocity u (assumed constant), the accel- 
eration of the fluid particle is given by 

a = -u x u x r - 2u x V (1) 


INTRODUCTION 

Recent trends in the design of high-speed machin- 
ery necessitated careful consideration of factors 
affecting gear durability. As a major failure mode of 
high-speed gearing, scoring is receiving significant 
attention. Scoring failures have been attributed to a 
sudden surge of the contact temperature of the gear. 
This contact temperature is directly related to the 
bulk temperature of the gear blank. Thus, an accurate 
prediction of scoring thresholds hinges on a correct 
evaluation of the bulk temperature. To evaluate the 
contact temperature, finite-element programs were 
developed by Patir (1977), Townsend (1981), and 
El-Bayoumy (1984) with the purpose of analyzing the 
steady-state thermal behavior of spur gears. The 
steady-state approximation required the assumption of 
constant heat-transfer coefficients. The numerical 
values of the coefficients were based on correlation 
with experimental measurements made by Townsend (1981). 
When operating conditions such as speed, applied load, 
initial temperature, and oil jet velocity, pressure, 
and viscosity were varied, correlation between theory 
and experiment was not always satisfactory. 

The difficulties arising from having to assume the 
proper values for these heat-transfer coefficients were 
overcome in El-Bayoumy (1985), where a transient ther- 
mal analysis approach was developed. The convection 
process in El-Bayoumy (1985) was based on an extended 
model from Van Heijningen (1974) that yielded a set of 
time-dependent coefficients. 

Closer examination of the hydrodynamic and thermal 
phenomena within the adjacent fluid film shows the film 
thickness remaining above the impingement point to 
decrease with time. The thermal boundary-layer thick- 
ness, on the other hand, grows steadily with time. An 
instant is therefore reached when the two thicknesses 
are equal. At that point the fluid film is fully satu- 
rated with heat and can no longer participate in the 
convection process. The wetted side of the gear tooth 
may therefore be divided into two regions: a convec- 

tive region, where the heat-transfer coefficients, 
derived from the extended Blok model, can be used and 
a nonconvective region with vanishing heat-transfer 
coefficients. 

The objective of this paper is to quantify the 
influence of the interaction between the thermal bound- 
ary and the oil film thickness on the transient thermal 
behavior of spur gears. In addition, the thermal wave 
propagation process of the hot spots are portrayed by 
instantaneous isotherms as they move through the gear 
tooth. The asymptotic approach towards steady state in 
the absence and presence of the thermal boundary layer 
is evaluated. 

Interaction Between Thermal and Hydrodynamic Effects 
Within the Fluid Film 

To evaluate how much of the adjacent lubricant 
actually contributes to the heat dissipation of the 
gear tooth, one needs to consider the hydrodynamic 
forces acting on a fluid particle near the tooth sur- 
face above and below the point of impingement .1 


1-The derivations given in this section are based on 
private communications with Dr. John Murdock of the 
Aerospace Corporation. 


where r is the radius vector to the particle, and 
V is the relative velocity of jet impingement. The 
two terms appearing on the right hand side of Eq. (1) 
are the centrifugal and coriolis acceleration compo- 
nents. Since r varies very little across the tooth, 
the centrifugal acceleration, for all practical pur- 
poses can be considered constant for points along the 


tooth profile. 


Its magnitude is 


V 2 /r 

8 


where 


V 

g 


is 


the velocity of the gear, at the pitch point the cen- 
trifugal acceleration as seen from the sign in Eq. (1) 
tends to strip off the oil layer from the surface of 
the tooth. The magnitude of the second term, the 
coriolis acceleration, is 2VgVj/r where Vj is the 
magnitude of the jet velocity relative to the gear. 
Since Vj is of the same order of magnitude as Vg, 
the coriolis component cannot be neglected. For 
points on the tooth profile below impingement, the 
direction of the coriolis acceleration is normal to 
the tooth surface pointing outward. Conversely, for 
points above impingement the coriolis acceleration 
points inward, as illustrated in Fig. 1, establishing 
a "pseudo-gravitational" field. Thus, the initial 
area A of convection heat flux q, characterized by 
the equation 


q = hA AT (2) 

covers the length of the tooth profile above the 
impingement point. How much of this area will partici- 
pate in heat dissipation can be determined when the 
thermal and fluid boundary layers are examined versus 
time. The thicknesses of these layers as a function 
of time are schematically illustrated in Fig. 2. The 
viscous layer develops in a time of the same order of 
magnitude as the oil spread time which is negligible 
compared with the heat-transfer time. While the vis- 
cous forces tend to attach the oil film to the tooth 
surface, centrifugal acceleration tries to strip it 
off. The amount of film remaining on the surface has 
been determined by Jefferis (1967), who showed that the 
thickness of this oil film is given by 



where v is the kinematic viscosity of the oil, x is 
the coordinate along the tooth profile measured from 
the impingement point (see Fig. 1), G is the centrifu- 
gal acceleration, and t is time. This relationship 
is depicted in Fig. 2 as a straight line with negative 
slope. 

The thermal boundary layer thickness is growing as 
a function of time according to the deWinter and Bloc 
(1974) formula 


H t = 1.362Vat (4) 

where a is the thermal diffusivity of the oil. The 
distance x^ from the point of impingement of the 
point at which the thermal boundary layer thicknesses 
of Eq. (3) is equal to the film thickness of Eq. (4) 
and is therefore given by 


1.85 — t 2 

V 


(5) 
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This is the point of heat saturation of the fluid film. 
For points on the tooth profile below this point (i.e., 
x £ x^), the fluid no longer contributes to the convec- 
tive dissipation of heat. In other words, point P]^ 
becomes the effective impingement point. For points 
above P^, a heat-transfer coefficient based on the 
solution of the one-dimensional heat equation 


3T 

3t 



( 6 ) 


is used. This coefficient may be written as 

hj = K/Uat) 1 / 2 (7) 

Mathematical Formulation of the Transient Thermal 
Problems in Gearing 

As noted in El-Bayoumy (1985), the partial differ- 
ential equation governing heat transfer within the gear 
sector shown in Figs. 3 and h is given by 


where 



C a 4- < It- 0 
3 y 2 a 3z 2 0 3t 


( 8 ) 





X = 


X 

R 
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and where 






kP , (T - T ) 

a o 


J 


(9) 


where 

n n/R 0 

n normal coordinate to the tooth profile 

hj convective heat transfer coefficient in tooth 

profile above the effective point of impinge- 
ment, W/hr-m 2 K (Btu/hr ft 2 *F) 

h s convective heat transfer coefficient on front 
and back of gear 

hj t convective heat transfer coefficient on top land 

h t convective heat transfer coefficient on remain- 

der of tooth including bottom land 

P^ dimentral pitch of gear, (in. - *) 

q* q/q a v 

Tj tooth involute boundary above saturation point 

T s front and back end faces of tooth 

T t remainder of tooth profile 

It was possible in El-Bayoumy (1985) to reduce the 
above boundary value problem through discretization to 
the recurrence relation 

(C/At + Hy)§ n+ 1 = [C/At - H(1 - 2)]§ n + P (11) 

where C and H are the capacitance and heat conduc- 
tivity matrices, 0 n is the dimensionless temperature 
vector at time t n , and 0 < y < 1. As noted in 
Zienkiewicz (1977), for numerical stability y > 1/2. 

Frictional Heat Vector Options 

The instantaneous heat generated during contact 
between the two gears may be written as 


T current gear temperature, °C (°F) 

T 0 initial gear temperature, *C ( a F) 

a gear diffusivity, m 2 /s (in. 2 /s) 

k thermal conductivity of gear material, W/m K 

(Btu/ft *F) 

R 0 outside radius of gear, m (ft) 

F face width of gear, m (ft) 

q av . average generated heat flux through the gear pro- 
file, W/hr (Btu/hr) 


q “ fe Wf | V s| (12) 

where 

J mechanical equivalent of heat, W/Nm (Btu/ftlb) 

W normal load intensity, N/m (lb/ft) 
b width of contact region, m (ft) 

f coefficient of friction 

V s instantaneous sliding velocity, m/s (ft/s) 


The associated boundary conditions are 


3© 

3n 




h . (s, t)R 

_J 9. 
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6 - P.R q* 
d o^ 


0 


on 



Assuming an even split between the pinion and the gear, 
the frictional heat input to either gear is 


q i 



(13) 
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I Various options from the literature (Anderson, 1981) 

} (10) are available for the calculation of the friction coef- 

ficient. The following choices are incorporated into 
the program (according to Kelly-Benedict, 1961): 


J 


0.0127 log 


3.17xl0 6 vA 

. FnV s V 2 J 


(UO 
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If the effect of surface finish is to be examined, a 
modified Kelly-Benedict formula is 



The results presented in this paper are based on the 
Kelley-Benedict formula (Eq. (1*0) as it correlates 
well with efficiency measurements given in Anderson 
(1981). 

RESULTS AND DISCUSSIONS 

The heat-transfer coefficients derived in Eq. (7) 
were used in El-Bayoumy (1985) for all points above the 
true impingement point. The computer program developed 
for the purpose of evaluating transient thermal behav- 
ior was mainly concerned with numerical stability of 
two- and three-point finite-difference schemes. The 
computer program was also developed for a comparative 
study of gear temperatures based on heat-transfer coef- 
ficients used in the steady-state program with those 
based on Eq. (7). Preliminary results were given in 
the form of snap shots of gear temperature contours at 
the initial stages of tooth engagement. Because the 
program does many iterations, it requires considerable 
time to reach the steady-state equilibrium. Because 
all the experimental data are at steady-state condi- 
tions, it was necessary to run the program for some 
time to match the experimental data. 

The case investigated here is that of a pair of 
eight-pitch standard spur gears with 28 teeth each run- 
ning at 10 000 rpm with a normal load per unit face 
width of 5912 N/cm (3378 lb/in.). The actual jet 
impingement depth is assumed at 87.5 percent of the 
tooth height. The resulting temperature history is 
displayed in three forms. First, at predefined time 
instants, a scan is made over all elements for the max- 
imum and minimum temperature rise. Then the range is 
divided into a user-defined number of equal-temperature 
increments. Each division is -thus associated with a 
temperature level. These levels can be presented in 
the form of contour plots such as those shown in 
Fig. 5. The analysis is first carried out neglecting 
boundary-layer interaction (assuming that the fluid 
film continues to absorb heat.) The jet is assumed to 
impinge 5* past the last point of tooth engagement. 

The temperature range at each instant is divided into 
10 equal divisions, or contour levels. Because the 
gear temperature will continuously vary with time, each 
contour level will define a different temperature level 
each time. For example, the first contour level defines 
a temperature rise of 26.4 ®C (47.5 °F) at 0.0167 msec 
(l a ), 35.1 °C (63.1 ®F ) at 0.05 msec (3°), and 8.5 *C 
(15.2 °F) at 0.167 msec (10°) after the beginning of 
tooth engagement. 

In order to better visualize the chronological 
development of the heat wave, contour levels are 
initially preassigned, and the progress of each contour 
is monitored step by step. Thus the same temperature 
history may be displayed in form of the fixed contour 
levels shown in Fig. 6. The range of the temperature 
rise during the first cycle of tooth engagement, having 
been determined as approximately 75 ®C (135 ®F), is 
divided into 10 equal divisions starting with the first 
contour level at 75 *C (135 ®F) and ending with 0 °C 
(0 ®F) for the last contour level. The progress of 
heat-wave propagation can easily be traced through the 
sequence of the snap shots as the gear tooth moves 
through engagement given in Fig. 6. First, lower level 
contours emerge then proceed through the gear tooth, 
followed by the appearance of higher level contours 
then a gradual withdrawal of these contours occurs as 
the contact nears the pitch point. The contours then 
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reemerge as the tooth mesh proceeds through the angle 
of recess until all the 10 contours are realized at the 
end of contact. A significant drop in temperature then 
follows due to the absence of frictional heat and the 
continued dissipation of heat through the fluid film. 

The results may be summarized in a single tempera- 
ture rise versus time graph (shown in Fig. 7) illus- 
trating how the maximum, average midplane, and average 
face temperatures on the loaded and coast sides of the 
tooth vary with time. As indicated in the figure, most 
of the thermal activity is limited to the duration of 
tooth engagement. Furthermore the highest temperature 
is highly localized at the tip of the driving gear and 
at a point between the lowest point of contact and the 
pitch point. In the beginning of tooth engagement the 
tip is significantly hotter than the other localized 
hot spot. 

The influence of boundary-layer interaction, as 
illustrated in Figs. 8 to 10, accentuates that differ- 
ence. This can be expected, because at the lower hot 
spot the fluid film continues to convect heat as the 
thermal boundary layer has not fully developed as it 
will when engagement approaches the tip. Thus, while 
the lower hot spot goes through a temperature rise of 
40.75 *C (73.35 ®F) during initial engagement, the tip 
can rise by as much as 96.02 °C (172.84 °F) (Figs. 8 
and 9) versus 77.24 °C (139.02 *F) (Figs. 5 and 6) if 
boundary layer interaction is neglected. Note that the 
lower hot spot is unaffected by boundary-layer interac- 
tion (compare Figs. 5(a) and 8(a)). 

When the simulation process is continued to steady 
state, significant differences in predicted tooth tem- 
peratures can arise when boundary-layer interaction is 
not neglected. This is clearly illustrated in Figs. 11 
to 16. The maximum temperature rise predicted with 
boundary-layer interaction is 119 °C (215 °F) (Figs. 14 
and 15) in close agreement with experimental findings 
shown in Fig. 12 from Townsend (1981), compared with 
78 °C (141 °F) (Figs. 11 and 13) without boundary-layer 
interaction. The errors are more significant when 
average tooth temperatures are compared. With boundary- 
layer interaction the average steady-state temperature 
rises 49 *C (88 *F) (Fig. 15). Without boundary-layer 
interaction a predicted average steady-state tempera- 
ture rise of only 14 *C (25 *F) is achieved (Fig. 13). 

In Fig. 16 the maximum temperature rise is plotted 
versus time over the entire simulation period. It is 
interesting to note that with boundary-layer effects 
the temperature attains a maximum before steady state 
is reached and scoring could take place even though the 
indicated steady-state temperature is below the scoring 
threshold of the tooth pair. 


SUMMARY OF RESULTS 

The interaction between the thermal and lubricant 
boundary layers has been incorporated into a transient 
thermal analysis program for predicting surface temper- 
atures in high-performance spur gears. As a result of 
this interaction a portion of the fluid film adjacent 
to the tooth becomes saturated with heat. This portion 
continues to spread with time until it covers the entire 
wetted area. Predicted temperature rise for a given set 
of gears was found to correlate well with previously 
measured data. If boundary layer interaction is 
ignored, the predicted temperature rise at the hottest 
spot could be off by as much as 42.5 percent. Two hot 
spots exist on the tooth surface of the driving gear: 
at the tip of the loaded side and below the pitch point. 
The gear reaches its hottest temperature before steady- 
state conditions are achieved. Through the use of pre- 
scribed temperature contours, the process of heat-wave 
propagation through the gear tooth can be animated. 
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FIGURE 1. - CORIOLIS ACCELERATION ABOVE AND BELOW IMPINGEMENT 
ON GEAR TOOTH. 
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FIGURE 2. - CHARACTERISTIC LENGTH AND TIME SCALES FOR OIL 
FILM ON ROTATING GEAR TOOTH. 
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FIGURE 3. - CONVECTIVE HEAT 
TRANSFER COEFFICIENT ZONES. 


FIGURE i|. - FINITE ELEMENT 
GRID. 
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FIGURE 5. - TEMPERATURE RISE CONTOURS (IN DEGREES CELSIUS) AT 
1° TO 46.5° OF GEAR ROTATION AFTER INITIAL TOOTH CONTACT. 
INITIAL NONEQUILIBRIUM CONDITIONS WITH 10 EQUAL CONTOUR LEVELS 
BETWEEN MINIMUM AND MAXIMUM TEMPERATURE. SPEED 10 000 RPM; 
LOAD 5912 N/CM (3378 LB/IN.); 87.5 PERCENT OIL JET IMPINGEMENT 
DEPTH; NO BOUNDARY LAYER EFFECT. 
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FIGURE 8. - TEMPERATURE RISE CONTOURS (IN DEGREES CELSIUS) AT 1° 
TO 48.5° OF GEAR ROTATION AFTER INITIAL TOOTH CONTACT. INITIAL 
NONEQUILIBRIUM CONDITIONS WITH 10 EQUAL PREASSIGNED FIXED CON- 
TOUR LEVELS BETWEEN MINIMUM AND MAXIMUM TEMPERATURE. SPEED 
10 000 RPM; LOAD 5912 N/CM (3378 LB/IN.); 87.5 PERCENT OIL JET 
IMPINGEMENT DEPTH; NO BOUNDARY LAYER EFFECT. 
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FIGURE 7. - TEMPERATURE RISE VERSUS TIME OR DEGREES ROTATION. INITIAL NONEQUILIBRIUM 
CONDITIONS. SPEED, 10 000 RPM; LOAD, 5912 N/cm (3378 LB/IN.); 87.5 PERCENT OIL JET 
IMPINGEMENT DEPTH; NO BOUNDARY LAYER EFFECTS. 
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FIGURE 8. - TEMPERATURE RISE CONTOURS (IN DEGREES CELSIUS) AT 1° 
TO 46.5° OF GEAR ROTATION AFTER INITIAL TOOTH CONTACT. INITIAL 
NONEQUILIBRIUM CONDITIONS WITH 10 EQUAL CONTOUR LEVELS BETWEEN 
MINIMUM AND MAXIMUM TEMPERATURE. SPEED, 10 000 RPM; LOAD 5912 
N/CM (3378 LB/IN.); 87.5 PERCENT OIL JET IMPINGEMENT DEPTH; 
WITH BOUNDARY LAYER EFFECT. 



FIGURE 9. - TEMPERATURE RISE CONTOURS (IN DEGREES CELSIUS) AT 1° 
TO 96.5° OF GEAR ROTATION AFTER INITIAL TOOTH CONTACT. INITIAL 
NONEQUILIBRIUM CONDITIONS WITH 10 EQUAL PREASSIGNED FIXED CON- 
TOUR LEVELS BETWEEN MINIMUM AND MAXIMUM TEMPERATURE. SPEED, 

10 000 RPM; LOAD, 5912 N/cm (3378 lb/in.); 87.5 PERCENT; OIL 
JET IMPINGEMENT DEPTH; WITH BOUNDARY LAYER EFFECTS. 
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FIGURE 10. - TEMPERATURE RISE VERSUS TIME AND ROTATION. INITIAL NONEQUILIBRIUM 
CONDITIONS. SPEED, 10 000 RPM; LOAD, 5912 N/cm (3378 lb/in.); 87.5 PERCENT OIL 
JET IMPINGEMENT DEPTH; WITH BOUNDARY LAYER EFFECTS. 
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FIGURE 11. - TEMPERATURE RISE CONTOURS (IN DEGREES CELSIUS) AT 1° 
TO 46.5° OF GEAR ROTATION AFTER INITIAL TOOTH CONTACT. WITH 10 
EQUAL PREASSIGNED FIXED CONTOUR LEVELS BETWEEN MINIMUM AND MAXI- 
MUM TEMPERATURE. SPEED, 10 000 RPM; LOAD, 5912 N/cm (3378 
LB/IN.); 87.5 PERCENT OIL JET IMPINGEMENT DEPTH; WITHOUT BOUND- 
ARY LAYER EFFECTS AFTER STEADY-STATE EQUILIBRIUM CONDITIONS. 
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FIGURE 12. - I.R. MICROSCOPE AND THERMOCOUPLE MEASURE- 
MENT OF GEAR TEMPERATURE VERSUS OIL JET PRESSURE FOR 
THREE LOADS, SPEED, 10 000 RPM; OIL JET DIAMETER, 
0.08 CM (0.032 in.); INLET OIL TEMPERATURE, 308 K 
(95 °F). 
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FIGURE 14. - TEMPERATURE RISE CONTOURS (IN DEGREES CELSIUS) AT 
1° TO 46.5° OF GEAR ROTATION AFTER INITIAL CONTACT AFTER 63 
GEAR REVOLUTIONS. WITH 10 EQUAL PREASSIGNED FIXED CONTOUR 
LEVELS BETWEEN MINIMUM AND MAXIMUM TEMPERATURE. SPEED, 10 000 
RPM; LOAD 5912 N/cm (3378 LB/IN.); 87.5 PERCENT OIL JET IMPINGE- 
MENT; WITH BOUNDARY LAYER EFFECTS. 
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GEAR ROTATION AFTER 63 REVOLUTIONS, deg 

FIGURE 13. - TEMPERATURE RISE VERSUS TIME OR GEAR ROTATION AFTER 63 REVOLUTIONS. 

SPEED, 10 000 RPM; LOAD, 5912 N/cm (3378 LB/ltO; 87.5 PERCENT OIL JET IMPINGE- 
MENT DEPTH: WITHOUT BOUNDARY LAYER EFFECTS AFTER STEADY-STATE EQUILIBRIUM CONDITIONS. 
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FIGURE 15. - TEMPERATURE RISE VERSUS TIME OR DEGREES ROTATION AFTER 96 REVOLUTIONS. SPEED, 
10 000 RPM; LOAD, 5912 N/cm (3378 LB/IN.); 87.5 PERCENT OIL IMPINGEMENT DEPTH AND WITH 
BOUNDARY LAYER EFFECTS AFTER STEADY-STATE EQUILIBRIUM. 
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FIGURE 16. - PROFILE MAXIMUM TEMPERATURE VERSUS GEAR REVOLUTIONS WITH AND WITHOUT BOUNDARY 
LAYER EFFECTS. 
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